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A versatile semiclassical approximation for intense laser-atom processes is presented. This uniform approxi- 
mation is no more complicated than the frequently-used multi-dimensional saddle-point approximation and far 
superior, since it applies for all energies, both close to as well as away from classical cutoffs. In the latter case, it 
reduces to the standard saddle-point approximation. The uniform approximation agrees accurately with numeri- 
cal evaluations for potentials, for which these are feasible, and constitutes a practicable method of calculation in 
general. The method is applied to the calculation of high-order above-threshold ionization spectra with various 
binding potentials: Coulomb, Yukawa, and shell potentials which may model Ceo molecules or clusters. The 
shell potentials generate high-order ATI spectra that are more structured and may feature an apparently higher 
cutoff. 



I. INTRODUCTION 

Sufficiently intense laser fields ionize atoms or molecules 
by the quantum-mechanical process of tunneling [[!]]. Both 
the tunneling process and the ensuing motion of the electron 
in the continuum are well accessible to semiclassical meth- 
ods. Tunneling generates a wave packet whose center follows 
a classical trajectory while the wave packet is spreading. It 
may or may not return to within the range of the ionic bind- 
ing potential. If it does, the well-known recollision-induced 
processes, such as high-order harmonic generation (HHG) or 
high-order above-threshold ionization (ATI), take place |[2| |3j]. 

In the tunneling regime, the quantum-mechanical transition 
amplitude can be analyzed, computed, and interpreted via the 
saddle-point approximation ft H]. Typically, the transition 
amplitude is represented by a multi-dimensional integral over 
the time t' at which the electron enters the continuum by tun- 
neling, the later time t at which it revisits the ion, and one 
or all components of the drift momentum k along its orbit 
in between those two times. For a specified final state, e.g., 
for given final momentum of the electron after the recollision 
or for given frequency of high-order harmonic emission, the 
saddle-point approximation selects those particular "quantum 
orbits" that contribute to this final state. These orbits are char- 
acterized by particular values of the parameters t, t' and k, 
which are complex numbers because of the tunneling nature 
of these orbits. For a specified final state, there are, in general, 
several contributing quantum orbits. Their contributions have 
to be added coherently, and this yields an interference pattern, 
which may appear very intricate, even though its physical ori- 
gin is simple 

Within the context of atoms in strong fields, the contribut- 
ing quantum orbits typically come in pairs. This may be best 
known from the Lewenstein model of HHG: For specified 
harmonic order within the "plateau", there are two quantum 



orbits whose contributions dominate the harmonic yield, the 
"long orbit" and the "short orbit". An electron on the long 
orbit starts earlier (by ionization) and returns later (for recom- 
bination) than an electron on the short orbit [Q]. This is a very 
general feature of intense-laser-atom processes and holds also 
for the more complicated orbits, which bypass the ion once 
or several times before the recombination process takes place 
[g|. For fixed laser intensity, the maximal HHG frequency 
or the maximal energy of an ATI electron obey classical lim- 
its [[| g], which are related to the maximal kinetic energy of 
the electron returning to the ion. For parameters approaching 
such classical limits, the two quantum orbits become more and 
more identical. If it were not for the fact that their parameters 
are complex, reflecting the birth of the electron by tunneling, 
the two orbits of a pair would coalesce at the classical cutoff 



The near coalescence of the orbits of a pair near a cutoff 
constitutes two problems for the saddle-point approximation: 
(i) treating the two saddle points as independent becomes an 
increasingly inaccurate approximation if they approach each 
other closely [l(J; should they actually merge into one, 
the standard saddle-point approximation diverges and hence 
is completely inapplicable; (ii) beyond the cutoff, in the clas- 
sically forbidden regime, both (complex) saddle points con- 
tinue to exist as formal solutions of the saddle-point condi- 
tions. One, however, has to be dropped from the transition 
amplitude, which is frequently, but not always, indicated by 
an exponentially exploding contribution. A rigorous analysis 
establishes that, actually, it is not possible to deform the con- 
tour of integration through such a saddle point by the method 
of steepest descent. In the framework of the theory of asymp- 
totic expansions, the global bifurcation of the steepest-descent 
contour from two visited saddles to a single visited saddle is 
known as the Stokes transition [11 |l2| ]. In previous work, 
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these problems were not treated in a systematic fashion. In 
this paper, we invoke a specific uniform approximation to 
solve both problems Jl3| ] . This turns out to be no more com- 
plicated than the standard procedure of treating the two saddle 
points as independent, because it uses exactly the same infor- 
mation as the standard procedure, namely, the values of the 
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action and its second derivative at the saddle points. Problem 
(i) is solved because the uniform approximation regularizes 
the saddle-point integrals close to the classical cutoff, while 
it reduces to the saddle-point approximation far away from 
the cutoff. Problem (ii) is solved by imposing the simple re- 
quirement of continuity on the transition amplitude, which 
automatically selects the appropriate branch of the multival- 
ued solution that does not contain the contribution of the un- 
physical saddle point beyond the classical cutoff. For a zero- 
range binding potential, the benefit of the saddle-point ap- 
proximation lies in the insight gained by the introduction of 
a few quantum orbits, which allow one to visualize the physi- 
cal mechanism behind recollision-induced processes. For the 
mere purpose of computation, the transition amplitude can be 
calculated as well, if not more easily, via a simple quadrature. 
We will use the zero-range potential as a test case, and find ex- 
cellent agreement for the uniform approximation, even where 
the usual saddle-point approximation fails. 

The zero-range potential is a valid model for the descrip- 
tion of a negatively charged ion in an intense laser field 
Jl4| , |l5[ [l6| ]. To what extent it can also be employed to 
model an atom in an intense laser field or, in other words, 
just how important the long range of the Coulomb potential 
is in this situation, has been the object of some debate. Sur- 
prisingly, it has turned out that at least for the qualitative ex- 
planation of most intense-field effects the Coulomb tail is not 
instrumental [Q, |5|, |l7p. Still more surprisingly, even the sub- 
tle quantum-mechanical enhancements of the ATI plateau at 
certain sharply defined intensities [Q are not specific to the 
Coulomb potential. In fact, a zero-range potential yields vir- 
tually the same enhancements, though at slightly different in- 
tensities [[9]]. 

From this point of view, being able to compare ATI spectra 
from zero-range and non-zero-range potentials is important. 
However, for non-zero-range potentials, a direct computation 
of the transition amplitude requires one to carry out a cumber- 
some multidimensional integral, and the uniform saddle-point 
approximation is the only viable approach. 

The purpose of this paper then is twofold. First, we deter- 
mine the specific uniform approximation that applies to the 
pairs of quantum orbits that appear in laser-induced rescatter- 
ing processes. Second, we use this uniform approximation to 
investigate the influence of the form of the binding potential 
on ATI. 

The plan of the paper is as follows: In Sec. II, we sum- 
marize the improved Keldysh approximation for the transition 
amplitude. In Sec. Q we discuss the saddle points that feature 
in the saddle-point approximation as well as in the uniform 
approximation, and review the saddle-point approximation as 
well as its problems close to classical cutoffs. In Sec. £v|, we 
determine the uniform approximation that overcomes these 
problems and describe its conceptual relation to the saddle- 
point approximation. In Sec. [V] we compare the ATI spec- 
tra obtained by these approximations to the numerical results 
for the zero-range binding potential. The uniform approxima- 
tion is then used in Sec. [Vl| to address the effect of a general 
(non-zero-range) binding potential on the ATI spectrum, us- 
ing Coulomb, Yukawa, and shell potentials as examples. A 



sum mary of the results and conclusions can be found in Sec. 
VII. 

We use atomic units (a.u.) throughout this paper. 



II. TRANSITION AMPLITUDE FOR RESCATTERING 
PROCESSES 

Strong-field phenomena, such as above-threshold ioniza- 
tion (ATI), are successfully described by transition amplitudes 
derived within a framework known as the strong-field approx- 
imation. This approximation neglects the binding potential in 
the propagation of the electron in the continuum, and the laser 
field when the electron is bound, which corresponds to treat- 
ing the process of rescattering in the first-order Born approx- 
imation on the background of the laser field. (The first-order 
Born approximation yields the exact differential cross section 
in the absence of the field both for the Coulomb potential as 
well as for the zero-range potential.) The ATI transition am- 
plitude for the direct electrons - electrons that leave the vicin- 
ity of the ion right after they have tunneled into the continuum 
- is the well-known Keldysh-Faisal-Reiss (KFR) amplitude 

/OO 
dt' (^(t')\v\Mt'))- (i) 
-oo 

The generalized transition amplitude, which includes one sin- 
gle act of rescattering, is given by $F\\\ 

/oo pt 
dt / df (^\t)\vu^ v \t,t')v\Mt'))- 
-oo J — oo 

(2) 

In both equations, V denotes the atomic binding potential, 
the final state is the Volkov state describing a charged parti- 
cle with asymptotic momentum p in the presence of a field 
with vector potential A(t), 

(r\^\t)) = exp (-'-J™ dr [p + A(r)] 2 ) e ^ +A ^ r , 

(3) 

and U( y,> (t,t') is the Volkov time-evolution operator, which 
describes the evolution of the electron in the presence of only 
the laser field. In Eq. ([[]), the electron, initially in the ground 
state W/jq (t')), is ionized into its final state at the time t 1 . In 
Eq. (0), an additional rescattering off the binding potential at 
the time t is accounted for. The amplitude (^J) incorporates 
the amplitude (jl|) for direct ionization in the limit where t' — > 
t. Hence, the two amplitudes must not be added [511]. The 
amplitude (Q) or closely related versions thereof have been 
used by several authors [|[ [lOj |2| §3j. 

If we insert the expansion of the Volkov propagator in terms 
of Volkov states, 

uW{t,H) = J d 3 k|4 v) (i)><v£V)|, (4) 
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into Eqs. (|l|) and (g), the transition amplitudes can be rewritten 

as 



dt' exp[iS p {t')] Vp , 



(5) 



and 



/ dt' / ^ke^^'-^FpuFko, (6) 

-oo «/— oo «/ 

where the corresponding actions are given by 

S p( t ') = -2j, dT [P + Mr)] 2 + \E \t (7) 



and 



1 f°° 

S p (t,t',k) = -- J dr [p + A(r)] 2 

- \j dT [k + A(r)} 2 + \E \t'. (8) 

The quantity \Eq | denotes the ionization potential of the atom. 
In this paper, we address the case of a linearly polarized 
monochromatic field, 



A(t) — Aoe x coscut, 



(9) 



with the ponderomotive energy Up = {A 2 (t)) t /2 = Aq/4. 

The representations (^|) and (^) are particularly useful if the 
form factors 

Vp k = (p + A(t)|V|k + A(t)) 

= / rf3r cx pH(p - k ) ' r M r ) (10) 

and 

Vuo = (k + A(i')|^|0) 

(-La / exp[-i(k + A(f )) • r]V(r)^(r) 



(11) 



can be calculated in analytical form. Within the strong-field 
approximation, the influence of the binding potential is en- 
tirely contained in these two matrix elements. For a zero- 
range potential, the form factors are constants. In this case, 
the five-dimensional integral (||) can be reduced to a one- 
dimensional integral over a series of Bessel functions, which 
can be readily computed numerically [ pl[ |24| ], In Sec. |v|, we 
will refer to the outcome of this procedure as the "exact re- 
sult". In general, however, a correspondingly "exact" evalua- 
tion of the matrix element (Q) has to deal with a multidimen- 
sional integral. 



III. SADDLE-POINT ANALYSIS 

For sufficiently high intensity of the laser field, correspond- 
ing to small Keldysh parameter 7 = y / \Eo\/2Up, ionization 



can be envisioned to proceed via the quasistatic process of 
tunneling [|5]| . The transition amplitudes (Q) and (g) are then 
conveniently computed via the method of steepest descent. 
Both the standard saddle-point approximation as well as the 
uniform approximation rest on this method, which approxi- 
mates the entire integral by the contributions from the vicinity 
of those points on the integration contour where the action 
is stationary, i.e., where the partial derivatives of the action 
with respect to the integration variables vanish. These points 
correspond to maxima of the integrand after a deformation of 
the original integration manifold, which is constructed such 
that the integrand decreases roughly like a Gaussian when one 
moves away from the vicinity of the saddles [pT|. 

In the current section, we first write down the equations that 
determine the saddle points, then describe the general proce- 
dure of identifying the relevant saddles, and finally discuss the 
saddle-point approximation. All these items are prerequisites 



for the discussion of the uniform approximation in Sec. IV 



A. Saddle-point equations 

For the rescattering amplitude (||), the saddle-point equa- 
tions are 



[k + A(i')f = 
[p + A(i)f = 



/ dr[k + A(r)] 



-2\Eo\, 
[k + A(t)f 

0. 



(12) 
(13) 

(14) 



Their solutions determine the ionization time t', the rescatter- 
ing time t, and the drift momentum k of the electronic orbit in 
between those two times, such that the electron acquires the 
asymptotic momentum p. Equations (|l2|) and (|l3]) are related 
to energy conservation at the ionization time and the rescat- 
tering time, respectively, and Eq. ([14]) determines the interme- 
diate electron momentum. For the direct amplitude (|j), only 
the ionization time t' need be determined, and the resulting 
equation is like Eq. ( |l2| ) with k replaced by the asymptotic 
momentum p. 

Evidently, Eq. ( |l2| ) has no real solutions t' as long as 
Eq 7^ 0, and in consequence t, t' and k are complex. Phys- 
ically, the fact that t 1 is complex means that ionization takes 
place through a tunneling process. The solutions (t, t 1 ) of the 
saddle-point equations for the linearly polarized monochro- 
matic field have been computed in Ref. They only 
depend on the ionization energy Eo and the photoelectron 
momentum p, but not on the shape of the binding potential, 
which enters the transition amplitude only via the form factors 
© and Q. 

A very important feature of the solutions is that they come 
in pairs. Let us denote the "travel time" by r = t — t'. Then, 
for given asymptotic momentum p and for the nth travel-time 
time interval nT/2 < Rer < (11+ l)T/2 (n = 1,2,...), 
there are two solutions having slightly different travel times. 
The parameters of two typical pairs of quantum orbits are dis- 
played in Fig. [lj 
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B. Classical cutoffs and Stokes transitions 

The original contour of integration in the amplitudes (||) or 
(^) is along the real axes, while the solutions of the saddle- 
point equations (p"2|)-(|l4}) are located off the real axes in the 
complex plane. A central question in the method of steepest 
descent then is, which of the various saddle points are vis- 
ited by the steepest-descent integration manifold. We shall 
call those the relevant saddle points. The steepest-descent 
manifold consists of pieces with a constant real part of the 
action. These pieces are glued together at zeros of the inte- 
grand, at which the phase of the action is not well defined. 
Usually, each piece visits only a single saddle point, which 
also determines the constant real part of the action. Only such 
pieces that are needed to connect the integration boundaries 
give contributions to the transition amplitude. The number of 
these pieces can change in a so-called Stokes transition, when 
two pieces merge at a certain value of a parameter (here we 
consider the photoelectron momentum p). On either side of 
the Stokes transition, the manifolds of the saddles of interest 
are glued together in different ways: on one side, both pieces 
are needed to connect the integration boundaries (plus, pos- 
sibly, other pieces related to different pairs of saddle points), 
while only one of the pieces is needed on the other. Note 
that in the latter case, too, there are still two solutions of the 
saddle-point equations, but only one of them is visited by the 
steepest-descent deformation of the original integration man- 
ifold [|<§. 

Merging of steepest-descent manifolds requires that the real 
parts of the actions of two quantum orbits become identical at 
a specific value of p, 



Re *Sp (ti , , k^ ) — R,e Sp(tj 1 tj 1 k 



(15) 



where i and j denote the saddle points of the given pair, and 
the times t s and t' s (s — i,j) depend on p. It follows from 
the physical mechanism behind high-order ATI that both sad- 
dles of each pair are relevant provided the asymptotic mo- 
mentum is classically accessible. For the pair of orbits hav- 
ing the shortest travel times (n = 1), this is the case if 
p 2 /2 < 10.007[/ P [0. The other pairs of orbits have smaller 
cutoff energies. 

The relevant saddle beyond the classical cutoff is the one 
that has the smaller imaginary part of the action at the Stokes 
transition []Z8|]. In the following we reserve the index i for 
this saddle. Saddle j only maintains a residual contribution to 
the transition amplitude after the Stokes transition, until it be- 
comes completely irrelevant in the so-called anti-Stokes tran- 
sition 



lmS p (ti,ti,ki) =!mS p (t j ,t' j , 



(16) 



The anti-Stokes transition coincides with the Stokes transition 
if both saddles actually coalesce. Otherwise, it frequently oc- 
curs very shortly after the Stokes transition. 

Exactly how the transition amplitude behaves close to the 
classical cutoff can only be described when the interplay of 
both saddles is taken into account in a systematic way, which 
is achieved by the uniform approximation. Before we turn to 



this approximation, we now discuss the standard saddle-point 
approximation. 



C. Saddle-point approximation 

Within the saddle-point approximation, the amplitudes (|5|) 
and (Q) are approximated by 



M. 



(SPA) 



E 



2m 



d 2 S p /dtj 



V p0 exp[iS p (t s )] (17) 



and 



= 5> s exp(^ s ), 

S 

S s — Sp (t s , t s , k s ) , 



A s = (2m) 



5/2 _ 



y p k„ Vk s o 



det S"(t,t',k)l 



(18a) 

(18b) 
(18c) 



respectively, where the index s runs over the relevant saddle 
points, and Sp(t, t' ,k)\ s is the five-dimensional matrix of the 
second derivatives of the action (|j) evaluated at the solutions 
of the saddle-point equations (|l2|)-([l4|). 

In explicit calculations, we will proceed slightly differently: 
First, we employ the saddle-point approximation to evaluate 
the three-dimensional integral over the intermediate momen- 
tum k in Eq. (Q), which enters the action (|8|) only quadrati- 
cally. This results in 



where 



dt 



— oo J — oo 



^V^'Vpk^^^o, (19) 



k(y,) = ~rb f t , dTA 



(r) 



(20) 



and S p (t,t') = S p (t, t' , k(i, t')). Then, we again make 
use of the saddle-point approximation to compute the two- 
dimensional integral over t and t' in Eq. (|r9|), which again re- 
sults in the amplitude (|l8]), where the actions and amplitudes 
are now computed by 



S s — S p (t s ,t s ), 



A. 



(2m) 



5/2_ 



Vn 



pk(t 3 ,t' 3 )Vk(t, ,t' s )0 

(t' s -t s )3 det S£(i,t') 



(21a) 
(21b) 



The corresponding saddle-point equations are Eqs. jl2l ) and 
([l3|) with k replaced by k(i, t'). Note that the values S s , A s of 
each saddle point are not changed, they are just obtained from 
a different set of relations in this more practical procedure. 

Upon approach to the classical cutoff, the two solutions that 
make up one pair come very close to each other. For an exam- 
ple, this is illustrated in Fig. [j]. The saddle-point approxima- 
tion (p"8|), however, treats different saddle points as indepen- 
dent. As mentioned in previous papers [||] and in the intro- 
duction, this leads to a quantitative and qualitative breakdown 
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of the standard saddle-point approximation near the cutoff of 
any pair of solutions, for two reasons: (i) This approximation 
can overestimate the contribution to the transition amplitude 
by several orders of magnitude (it actually diverges if both 
saddles coalesce), (ii) In previous papers, the spurious sad- 
dle has been dropped after the classical cutoff by requiring a 
minimal discontinuity of the transition amplitude. Still, the 
discontinuity remains finite and noticeable. 

A smooth suppression of the spurious saddle can be 
achieved if both quantum orbits are well separated at the 
Stokes transition (which is, however, not the case for physi- 
cally accessible parameters in ATI), by a regularization that 
has been derived in the general framework of asymptotic ex- 
pansions jl^]. Thereby, the contribution of the spurious sad- 
dle is suppressed by multiplication with the error function 

erfc(-i/) = — / drexp(-r 2 ), (22) 

with the argument given by 

v = , (23) 

^2|Im[S p (^)-S p (^)]| 

The argument v vanishes at the Stokes transition ( |l5| ) and di- 
verges at the anti-Stokes transition (|l(|), after which the spuri- 
ous saddle drops out completely. Note that this automatically 
prevents an exponential growth of the amplitude of the spu- 
rious saddle in the approximation (|8j), because the saddle is 
dropped while the imaginary part of the action is still positive 
(namely, equal to the imaginary part of a physical saddle). 

This regularization procedure is not accurate enough in 
the present problem because the Stokes transitions take place 
while the saddles are not sufficiently separated (cf. Sec. |y|). 
On the other hand, the Stokes transitions are already built into 
the uniform approximation, to which we turn now. 

IV. THE UNIFORM APPROXIMATION 

The saddle-point approximations ( |l7| ) and ( |l8| ) are obtained 
by expanding the action function S p to second order in the 
integration variables about each saddle point, and then do- 
ing the ensuing Gaussian integrals. These approximations are 
valid if the expansion of the action holds until the integrand 
has become much smaller than it was at the saddle point, so 
that the integration can be extended to infinity. The saddle- 
point approximation breaks down when the difference of ac- 
tions | Si — Sj | of two quantum orbits with similar coordinates 
becomes of order unity, such that the expansion about saddle 
point i becomes inaccurate close to the saddle point j, and 
vice versa. For the quantum orbits in ATI this happens when 
the energy approaches the classical cutoff. The remedy of- 
fered by the theory of asymptotic expansions is to improve 
the expansion of the action function in the neighborhood of 
saddles i and j by including higher orders in the coordinate 
dependence and to take the resulting approximate integral as 
a collective contribution of both saddle points. 



What is often not observed is that the resulting uniform ap- 
proximation can be written in such a form that no additional 
information on the quantum orbits is needed, i.e., the cumber- 
some expansion in the coordinate dependence actually can be 
circumvented. The derivation proceeds in two steps. First, we 
write down the so-called diffraction integral which describes a 
pair of orbits which might be close to each other or well sepa- 
rated. Then, we determine the parameters of the formal expan- 
sion in terms of the quantities that enter the standard saddle- 
point approximation, from the observation that the conven- 
tional saddle-point approximation ( p"8| ) has to be recovered in 
the limit where the saddle points are sufficiently well sepa- 
rated. 

For the first step, we observe that it is precisely two quan- 
tum orbits that closely approach each other near each cutoff. 
According to the splitting lemma of catastrophe theory [p^], 
the parametrization of the integration domain can be rectified 
such that the orbits approach each other along one of the (ap- 
propriately chosen) coordinate axes (denoted by x in the fol- 
lowing). This is the only direction where higher orders in the 
coordinate expansion of the action have to be included, while 
the expansion in the other coordinates can be restricted to sec- 
ond order such that these can be integrated out by the usual 
saddle point approximation [this is similar to integrating out 
k in the transition from Eq. ( |l8| ) to Eq. (|T])]. Hence the con- 
tribution of the pair of quantum orbits (denoted by i and j) 
to the transition amplitude can be reduced, in principle, to a 
one-dimensional diffraction integral of the general form 

M i+j = j dxg(x)exp[iS(x)], (24) 

where the action accounts for these two saddle points and 
the integration boundaries c u , c\ in (complex) infinity are as- 
sumed such that the integrand decays to zero and the inte- 
gral converges. Moreover, an expression that reduces to the 
conventional saddle-point approximation when the quantum 
orbits are well separated will be obtained if we allow for a 
linear coordinate dependence in the function g(x). This moti- 
vates the use of the normal forms (for a derivation in another 
semiclassical context, see Ref. [[l3||) 

S(x) = S + ex — ax 3 , g(x)=g n +gix. (25) 

Here we have chosen the origin of the coordinate system ex- 
actly in the middle between the two saddles, which have co- 
ordinates Xij = ±y£/3a and coalesce when e = 0. 

The uniform approximation that we introduce here differs 
from an earlier regularization method [pt [Tcfa, where the action 
was expanded to cubic order about the stationary point corre- 
sponding to the classical cutoff. This led to the absence of the 
linear term in the function g(x) in Eq. (]25j). It is precisely this 
term whose presence allows us to match the standard saddle- 
point approximation both near the cutoff and away from it. 
Thus, the method of Ref. [[)]] coincides with the uniform ap- 
proximation near the classical boundary, but deviates from it 
and from the exact solution farther away from the cutoff re- 
gion. 
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FIG. 1: Saddle points as a function of energy for a Keldysh parameter of 7 = 0.975 and scattering angle 9 — 0. The first, second and third 
column give the start time, the return time, and the intermediate drift momentum, respectively. The panels present the paths in the complex 
plane that are followed by the saddle points as a function of the final energy, which is indicated by the numbers, which are in multiples of Up. 
The upper row gives the saddle points for the pair of orbits with the shortest two travel times (1 + 2), the lower row for (9 + 10), which is one 
of the pairs with the longest times considered in this paper. The figure shows how the saddle points of a pair approach each other most closely 
near the classical cutoff. In each case, the contribution of the orbit that is drawn dashed is dropped after the cutoff. 



With expansion (|25|) inserted into the original integral 
the amplitude Mi + j reduces to a sum of Bessel functions, 



M i+ j = ^2-kAS/Z exp(iS + wr/4) 

x {A[J 1/3 (A5) + J_ 1/3 (AS)] 
+ AA[J 2/3 (AS) - J_ 2/3 (AS)]}, 
AS = (Si-Sj)/2, § = (S i + Sj)/2 1 
AA = {A l ~iA J )/2, A = (iAi - Aj)/2, 



(26) 



where the four independent parameters S, AS — 
2e 3 / 2 (27a)" 1 / 2 , A = 5o (-2«) 1 / 2 a 1 / 4 (3e)~ 1 / 4 , and A^ = 
gi(27Tz) 1 / 2 e 1 / 4 (3a)~ 3 / 4 have been expressed by the ampli- 
tudes and actions that result from the saddle-point approxima- 
tion of the diffraction integral ([24]). 

The uniform approximation is defined by inserting into Eq. 
(^) the actions and amplitudes < Jl~8| ) of the respective pair of 
quantum orbits (which we denoted by i and j). We wish to 
stress that it is not necessary to obtain the expansion param- 
eters S, e, a, g , and gi by explicitly carrying out the expan- 
sion (|2^). Indeed, knowledge of the explicit dependence on 
these parameters is not even desired because it can be ma- 



nipulated by a coordinate transformation, while the original 
integral is invariant under smooth changes of the coordinate 
system. For the saddle-point approximation (JT^), invariance 
with respect to coordinate transformations is ensured trivially 
for the actions S s , while the amplitudes A s are invariant be- 
cause the Jacobian of a transformation contributes a factor to 
g which is cancelled by t he d etermi nant of second derivatives 
of the action, see Eqs. (18c) and (2 It). This is the reason 
why we express the expansion coefficients in Eq. ( ]25| ) by the 
coordinate-transformation invariant quantities Sij of the 
saddle points. Indeed, it is a simple exercise to verify with the 
help of the asymptotic behavior 



J±u{z) ~ — 
7TZ 



1/2 



cos(z =F vtt/2 — 7r/4) 



(27) 



of the Bessel functions for large z that the saddle-point ap- 
proximation (|l8]) is recovered from the uniform approxima- 
tion (H) in the limit of large AS. 

Finally, let us demonstrate that the uniform approximation 
is also capable of describing the Stokes transition, in which 
one of the two saddles is rendered irrelevant. The Bessel func- 
tions in (|26|) assume complex arguments and are multi-valued 



7 



functions, depending on the integration contour taken in their 
integral representation. The functional branches can be dis- 
tinguished by the number of saddles which are visited by a 
steepest-descent deformation of the contour, in complete anal- 
ogy with the procedure for the original integral (ph. Hence, 
when the condition (jlj) is fulfilled one not only observes a 
Stokes transition in the original integral, but also encounters 
a Stokes transition in the defining integral of the Bessel func- 
tions. The proper branch of the function will automatically 
be selected by requiring a smooth functional behavior. The 
choice of branches beyond the Stokes transition corresponds 
to replacing the Bessel J functions by Bessel K functions, 



M, 



i+3 



^/2iAS/iTexp(iS) 
x [AK l/3 (-iAS) 



■iAAK 2/3 (-iAS)] . (28) 



From the usual asymptotics 



/ 7T \ V 2 



(29) 



of the Bessel K function for large z one verifies that in this 
case only saddle i contributes to the saddle-point approxima- 
tion. 

In summary, in the uniform approximation the sum of 
saddle-point amplitudes dlif ) of each pair of quantum orbits 
is simply replaced by the collective amplitude (p6|). The uni- 
form approximation improves the saddle-point approximation 
such that it works even when two quantum orbits approach 
each other so closely that one cannot locally expand about ei- 
ther one, as is the case close to their classical cutoff. It also 
works well far away from classical cutoffs, because it includes 
the saddle-point approximation as a special case which is re- 
covered for \AS\ > 1. This can happen in two ways: (i) 
when the saddle points become well separated as a system pa- 
rameter (such as p) is varied, or (ii) in the strict semiclassical 
limit when for fixed system parameter the Keldysh parameter 
is decreased (given AS ^ 0). Also, the Stokes transition at 
the classical cutoff is automatically built into the uniform ap- 
proximation. Most notably, the uniform approximation is of 
the same practical simplicity as the saddle-point approxima- 
tion since it involves the same amplitudes A s and actions S s 
defined in Eqs. (p"i|). 



V. COMPARING THE VARIOUS APPROXIMATIONS 

In this section, for the zero-range potential we compare the 
approximations discussed in the previous sections with the ex- 
act integration of Eq. (||). First, let us consider ATI spectra in 
the direction of the electric field of the laser. Such a spectrum 
is composed of the contributions of direct and of rescattered 
electrons. The former quickly decrease after their classical 
cutoff at 2Up. The latter form an extended plateau with its 
classical cutoff at 10[/p, whose yield is below that of the di- 
rect electrons by several orders of magnitude. The cutoff at 
10Up is related to the pair of orbits with the shortest travel 
times. The other pairs of trajectories, which have longer travel 
times, have cutoff energies below this value (see, e.g., Ref. [0] 
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FIG. 2: Photoelectron spectra for a zero-range binding potential and 
Up/co — 3.58, W = 0.073 a.u., and a ground-state energy of 
Eo = —0.5 a.u. The spectrum is in the direction of the electric 
field of the laser, 8 — 0. Part (a) shows spectra computed using the 
saddle-point and uniform approximations, compared with the pho- 
toelectron yield obtained by computing the integral ([j) exactly. We 
take into account the two direct trajectories and five pairs of rescat- 
tered trajectories. The approximate energy positions of the Stokes 
transitions, which coincide with the respective classical cutoffs, are 
indicated by arrows. Part (b) displays spectra computed by means of 
the uniform approximation, for direct, rescattered, and both types of 
electrons, and compares these with the exact integration. 



for a more complete discussion). In the figures that follow, we 
consider up to 5 pairs of electron trajectories, those with the 
shortest travel times. To each trajectory, we associate a pos- 
itive integer number which increases with the corresponding 
travel time. 

The outcome of this comparison is displayed in Fig. ^a). 
In general, there is a good qualitative agreement between the 
saddle-point approximation and the exact solution (note, how- 
ever, that the scale is logarithmic in this figure.) Quantita- 
tively, however, there are marked discrepancies, which occur 
in those energy regions where the saddle points that consti- 
tute a particular pair approach each other and can no longer 
be treated as independent. 

In previous work [0], the unphysical contribution of one 
of the saddle points was eliminated by hand as soon as the 
energy crossed the Stokes line (|l5|). This causes the cusps 
in the spectra, which can also be seen in Fig. ^|(a). This is 
not very satisfactory, since the discrepancies in the ATI signal 
may amount to almost one order of magnitude. This problem 
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caused by the trajectories with still longer travel times that 
have not been included. Those do not contribute for larger an- 




gles. The saddle-point approximation, on the other hand, ex- 
hibits large discrepancies with the exact results near the clas- 
sical cutoffs. For the chosen photoelectron energy of 8.01C/p, 
there are only three relevant cutoffs, corresponding to the pairs 
of trajectories 1+2, 5+6, and 9+10. The remaining pairs of tra- 
jectories do not contribute, since their cutoffs are significantly 
below8.01[/ P . 



e (degrees) VI. INFLUENCE OF THE POTENTIAL ON 

RESCATTERING PROCESSES 



FIG. 3: Angular distributions of photoelectrons for the zero-range 
potential case, computed with the saddle-point and uniform approx- 
imations, compared to the exact yield. The field parameters are 
Up/co — 35.8, u) = 0.0584 a.u., the ground-state energy is chosen 
as Eq — —0.9 a.u., and the photoelectron energy is e — 8.01f/p. 
The angles of Stokes transitions are marked with arrows. 



is particularly critical if the intensity of the driving field is not 
so high. In this case, the various cutoff energies are relatively 
close to each other, so that the artifacts affect a broad energy 
region. Thus, a more accurate approximation is desirable and 
even necessary, in case the integral (Q) cannot be carried out 
exactly, as is the case for any potential other than the zero- 
range potential. 

One possibility to eliminate such effects, shown in Fig. ||(a), 
is the Stokes regularization, Eq. (p2|). This smoothes out the 
cusps, without, however, eliminating them completely. 

Far superior results are obtained by the uniform approxi- 
mation, given by Eqs. ( p6| ) and (p8|), respectively. The spec- 
trum computed in this way almost perfectly agrees with the 
exact result. The remaining differences between the uniform 
approximation and the exact integration occur near the inter- 
ference minima and are due to the contributions of pairs of tra- 
jectories with longer travel times that have not been included. 
This is indicated by the minor differences in the spectra com- 
puted with the uniform approximation using 3 and 5 pairs of 
trajectories, cf. Fig. ^b). 

Figure ||(b) shows that the exact spectrum is well repro- 
duced by the uniform approximation for all energies. The 
figure also separately displays the contribution of the direct 
electrons p0|]. One observes that interference between the 
rescattered and direct electron trajectories is only important 
within a small energy region, between AUp and 6Up pT|]. 
Above and below this energy range, either the rescattered or 
the direct electrons completely dominate the spectrum, so that 
interference only leads to minor effects. 

The superiority of the uniform approximation over the 
saddle-point approximation becomes particularly impressive 
if spectra are displayed on a linear scale. This is done in 
Fig. U for an angular distribution at fixed energy. Both with 
the saddle-point approximation and the uniform approxima- 
tion, the 10 shortest trajectories are considered. The uniform 
approximation, again, yields excellent agreement with the ex- 
act result. Minor differences, for small scattering angles, are 



The preceding section has shown that the uniform approx- 
imation is a very dependable method, yielding results very 
close to those obtained from the exact integration. The lat- 
ter, however, is only feasible for a binding potential of zero 
range. Therefore, we will rely on the uniform approximation 
to investigate how the form of the binding potential affects 
the photoelectron spectrum. The transition amplitude (||) was 
derived in the context of one electron bound by the potential 
V(r). In order to simulate a many-electron atom, it can be 
reasonable to use in the transition amplitude (Q) different po- 
tentials V(r) for the electron when it tunnels out and when it 
rescatters [ 10 1 . In Refs. pC| , |23p, the effect of the rescattering 
potential on the general shape of the high-order spectrum and 
the ratio of direct over rescattered electrons were investigated 
as a function of the applied field, for the pair of the two short- 
est orbits. In particular, the dependence on the atomic species 
was modeled by a Thomas-Fermi potential. Here, for various 
model potentials, making use of the additional power afforded 
by the uniform approximation, we will concentrate on the de- 
tailed shape of the angular-resolved energy spectrum and on 
the contributions of the orbits with longer travel times. 

Throughout, we shall use the results for the zero-range po- 
tential 



V(r) = _gL=5(r)|-r 



(30) 



y / 2\EV\ 'dr 
as a benchmark. Its form factors 

© and @ are constants 



Vi- 



and 



P k 



Vko 



(2^) 2 V / 2^oT 
(2|£o|) 1/4 



2tt 



(31) 



(32) 



A. Influence of the Coulomb tail 

In this subsection, we investigate the influence of the long- 
range Coulomb potential on above-threshold ionization. This 
is particularly interesting since for hydrogen ATI spectra have 
been extracted from a high-precision numerical solution of the 
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time-dependent Schrodinger equation (TDSE) [|32|], so that we 
can compare the strong-field approximation with an exact so- 
lution. 

The form factors of the Yukawa potential V(r) = 
— Zexp(— ar)/r are 



Vr 



pk 



2tt 2 (p - k) 2 



(33) 



and 



Vko = 



V2 



Z b ' 2 



7T (Z + a) 2 + [k + A(i')] 2 
V2 Z 5 / 2 



7T (Z + a) 2 -2|£ | 



(34) 



where the saddle-point equation ( |l2| ) has been used in the last 
line. Hence, in the saddle-point approximation, Vko acts as 
a constant; indeed, this is the case for any spherically sym- 
metric potential. This constant determines the total ionization 
rate, but has no effect on the shape of the spectrum. Another 
consequence is that the spectrum of the direct electrons, de- 
scribed by the amplitude (Q), is independent of the form of the 
binding potential because it only depends on VkOi m contrast 
to the spectrum of the rescattered electrons. 

The Coulomb form factors can be retrieved from Eqs. (B3b 
and ( J34] ) in the limit a — > 0. Since in this case E = —Zy2, 
this leads to the well-known divergence of the Coulomb form 
factor ( p4| ) [0]. This has no effect on the shape of the spec- 
trum, and the absolute scale can be reestablished, too [p3|. 

In Fig. [|, we compare ATI spectra for the zero-range, the 
Yukawa, and the Coulomb potential. In view of the Coulomb 
divergence of Vko we used the zero-range form factor ([32]) for 
all potentials |3~4j|. As expected from Eq. ([33]), there is a sup- 
pression of the photoelectron yield for the higher energies in 
the Coulomb and Yukawa cases. This effect is present for all 
pairs of trajectories. For the Coulomb potential, there is an 
additional enhancement of the rescattered yield for low ener- 
gies, which does not occur in the zero-range or short-range 
cases. This enhancement is due to the functional form of Vpk- 
Clearly, if the screening parameter is small enough, this ef- 
fect is also present for the Yukawa potential. Furthermore, for 
these latter potentials, there is a reduction in the plateau in- 
tensity as the screening parameter is increased. Evidently, the 
form factor ( |33| ) for the Coulomb potential always exceeds the 
one for the Yukawa potential. 

The parameters of Figs. |2| and | correspond to those cho- 
sen in Ref. p2]], where the results of a numerical solution 
of the three-dimensional time-dependent Schrodinger equa- 
tion for hydrogen are reported and ATI spectra are extracted 
from the former. The agreement between the Coulomb result 
of Fig. [|(a) and Fig. 2 of Ref. is good and even quantita- 
tive. We notice that the pronounced dip in the spectrum near 
8Up, which is due to destructive interference of the contribu- 
tions of the shortest two orbits [cf. Fig. ^ (b)], is almost at 
the same position in both calculations. The next destructive- 
interference minimum from these two orbits occurs just below 
6Up. The contributions of the longer orbits [cf. Fig. [j](c) and 




Photoelectron energy (units of U p ) 



FIG. 4: Photoelectron spectra for the zero-range potential, compared 
with those for the Coulomb and Yukawa potential, for the same field 
and atomic parameters as in Fig. [2] Panel (a) shows total spectra, 
while panels (b) to (d) exhibit the contributions of individual pairs of 
rescattered orbits. 



(d)] partially fill in this minimum, leaving only a shoulder in 
the total spectrum (a). The exact calculation [ p2[ ] features a 
slightly more pronounced minimum at the same position. Re- 
markably, the two interference minima in the total spectrum 
at low energy near 0.5Up and 2Up, which are due to the di- 
rect electrons and the amplitude (||), are also clearly reflected 
in the exact calculation [[32p at about the same positions. The 
overall drop of the spectrum from the direct electrons to the 
final maximum of the rescattered electrons preceding the cut- 
off is more pronounced in the exact calculation by about half 
an order of magnitude [ f3~5[ ] . 

In Fig. ||, we investigate the ATI spectra for several screen- 
ing parameters a of the Yukawa potential. In this figure, we 
also address the question of how the form factor Vko affects 
the photoelectron yield. The figure clearly shows a global 
shift in the photoelectron signal, which increases for decreas- 
ing a. In this sense, our results are in agreement with those in 
Ref. [^2j. It is, however, not expected that this yield increases 
indefinitely. In fact, its limit for a — > should be given by the 
TDSE results p2|]. Because of the singularity for hydrogen 
in Vko for vanishing screening parameters, such a comparison 
is beyond the scope of the strong-field approximation. Addi- 
tionally, there is an enhancement of the photoelectron yield 
at lower energies, similar to those occurring in the Coulomb 
case, which disappears as a is increased, which is in agree- 
ment with the previous figure. 
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FIG. 5: Photoelectron spectra for the Yukawa potential, the same 
field and atomic parameters as in the previous figure, and several 
screening parameters a. Part (a) shows the resulting spectra for 
the direct electrons and the five shortest pairs of rescattered orbits, 
whereas part (b) shows the contributions from the shortest pair of 
rescattered trajectories. 



B. Shell potentials 

Spherical shell potentials have been used for modeling clus- 
ters or molecules such as C§q. Recently, ATI has been ob- 
served experimentally for Ceo in the direct-electron energy re- 
gion [[36]]. Therefore, in this section we investigate how such 
potentials affect the ATI spectra in the direct and rescattered 
regions. Let us first consider a spherical (5-shell, 



V(r) =-V 6(r-r ), 



with 



V 



y/2\Eo\ 



1 - exp[-2 v /2^r ] ' 



(35) 



(36) 



where Eq again denotes the binding energy of the ground 
state. Ionization from such a potential was investigated in the 
past [p7], for weaker laser fields. The corresponding form 
factors (10|) and ([H]) are 



Vp k = — 



V r 



: sin[^/ (p - k) 2 r ] 



(37) 



FIG. 6: Photoelectron spectra for the shell potential (p5[), compared 
with the zero-range case. The ionization potential was taken as 
\Eo\ = 0.274 a.u. and the cluster radius as ro = 6.7 a.u. The 
field parameters are la = 6.5 x 10 13 W/cm 2 , and u = 0.057 a.u. 
This yields an excursion amplitude of ao = 13.2 a.u. and a Keldysh 
parameter 7 = 0.9805. In part (a) we take into account the five short- 
est pairs of trajectories, whereas in part (b) only the shortest pair is 
considered. 



and 



V k0 



V C 



respectively, with 



C 



exp(2^/2\EV\r ) - 1 - 2^2\Eo\r 



1/2 



(38) 



(39) 



For the (5-shell potential, Vpk is an oscillating function, and 
Vko is a constant as always. Thus, in the following, we con- 
centrate on the influence of V p k on the resulting spectra. We 
consider typical Ceo parameters, taken from Ref. |36|. The 
external field is chosen such that its intensity is still below the 
Ceo fragmentation threshold, but the electron excursion am- 
plitude [ p8| ] is roughly twice as large as ro. Furthermore, the 
Keldysh parameter is about unity. Thus, the rescattering pic- 
ture is still expected to be applicable. 

In Fig. [| we compare the photoelectron spectrum for the 
(5-shell and for the zero-range potential, within the uniform 
approximation. In order to assess the efficiency of rescatter- 
ing, in either case we used for Vko the zero-range result (|32|). 
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The figure shows that the (5-shell potential rescatters more ef- 
ficiently than the zero-range potential by about one order of 
magnitude. If the form factor ( |38| ) is taken into account, an 
additional global increase in the yield occurs. However, in 
the (5-shell case, the rescattering plateau on the average has a 
downward slope, in contrast to the zero-range case where the 
slope goes up. 

The most interesting feature, however, is that the rescattered 
spectrum of the (5-shell potential is much more structured than 
it is for the zero-range potential, with several additional oscil- 
lations. Such oscillations are due to the form factor (|37|), and 
are already present for the contributions of the shortest pair of 
trajectories, as shown in Fig. [Jb). An unexpected side effect 
of these oscillations is the effective increase of the plateau cut- 
off energy by about two units of Up for the shell versus the 
zero-range potential, which can be observed in Fig. [| Since 
the laser intensity is the same in both cases, the rescattering 
cutoff would be expected at the same energy, too. However, 
the shell form factor has a zero around the energy of 9.5Up, 
where the zero-range spectrum features its final maximum. 
This moves the final maximum of the shell-potential spectrum 
up to a higher energy. 

In order to investigate these oscillations in more detail, in 
the following we will look at contributions of individual tra- 
jectories to the photoelectron yield for the (5-shell, in compari- 
son to the zero-range potential. Since the uniform approxima- 
tion requires pairs of trajectories, we will use the saddle-point 
approximation for that purpose. Whenever dealing with a pair 
of trajectories, we will consider the uniform approximation. 

Figure ^ displays these results, for several rescattered tra- 
jectories. In case V p k is constant, as is the case for the zero- 
range potential, all oscillations present in the spectra come 
from interference terms. The contributions of individual tra- 
jectories are nearly constant in the classically allowed regime 
and do not produce any substructure. For the (5-shell, however, 
V^,k is oscillatory and produces its own maxima and minima 
in the spectrum. However, comparing Fig. ^a) and (b) we 
observe that the contributions of the longer orbits tend to re- 
store the minima of the shell-potential spectrum to those of the 
zero-range. Only the highest-energy minimum near 9.5Up is 
left unaffected, since the longer orbits do not contribute to this 
energy. 

In particular, the minima are given by Re yjj» — k) 2 = 
nir/ro, where n is an integer. To a first approximation, the 
drift momentum k can be neglected with respect to the mo- 
mentum p, so that the energy positions of the minima, in units 
of the ponderomotive energy, are roughly given by 



P 



2U P r z U P 



(40) 



This expression is expected to work better for longer excur- 
sion times, since, accordingto Eq. (|I4|), k oc l/(t — t') . This 
can already be seen in Fig. |J where the saddle points as func- 
tions of the energy are depicted. For a pair of trajectories with 
short travel times, the start and the return times, as well as the 
intermediate momentum k, vary considerably with the pho- 
toelectron energy. For a long travel time, on the other hand, 
these quantities are nearly constant, in the classically allowed 
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FIG. 7: Contribution from individual trajectories to the rescattered 
photoelectron spectrum for the shell potential, in comparison to the 
zero-range case. We consider the same parameters as in the previous 
figure. The labels i and j refer, in part (a), to the third-shortest pair, 
denoted by (5+6), and in part (b) to the eighth-shortest pair, denoted 
by (15+16). For the terms \i + j\ 2 we applied the uniform approxi- 
mation. For the terms |i| 2 + |j| 2 , we applied the Stokes regularization 
jjj) to the diverging trajectory. The dashed vertical lines in the fig- 
ure separate the classically allowed and forbidden energy regions for 
the respective orbits. The dotted and dashed gray lines in part (a) 
denote the individual contributions of 5 and 6, respectively, for the 
zero-range case. 



region. Furthermore, the return time, as well as the intermedi- 
ate momentum, are almost real and k is very small. 

Clearly, there exist deviations from Eq. (E7j) due to the fact 
that k is non-vanishing and complex, t and rare complex, and 
due to the time dependence of the intermediate momentum. 
For instance, a feature that is not explained by Eq. ( f+(i| ) is a 
shift in the oscillations of the longer trajectory, with respect 
to those of the short one. This feature occurs for all pairs 
of trajectories, and decreases as the travel times get longer. 
A qualitative estimate of these deviations can be obtained by 
considering y/(p — k) 2 up to first order in %/k 2 , and the pair 
and (t2,f 2 ) — (ti — e,t\ + e') up to first order in 
e, e'. This gives a shift in the minima, which is proportional to 
s/(t — t'), confirming the results presented in Fig. f^. 

Now we turn to other shell potentials. Similar results are 
obtained for a more realistic square well, of the form V(r) — 
—Vo for n < r < r 2 , and zero otherwise. Since, in nature, the 
sharp edges present for a (5-shell or a square well are smoothed 
out, it is of interest to investigate whether the additional oscil- 



12 



2 





o 

CL 

O 

8> 




2 4 6 8 10 12 
Photoelectron energy (units of U ) 



FIG. 8: Contribution from the shortest pair of trajectories to the 
photoelectron spectrum for the Gaussian potential (J4- 1|>, compared to 
the 5-shell case, for several widths a and the same parameters as 
in the previous figure. Parts (a) and (b) depict |1 + 2| 2 and |2| 2 , 
respectively. The prefactor Vo, for \Eo\ — 0.274, was computed by 
solving the time-independent Schrodinger equation numerically. 



lations are also present for smooth potentials that approximate 
Eq. (35). One such example is the Gaussian potential 



V(r) =-V cxp[-(r-ro) 2 /<T 2 



(41) 



For vanishing width, we recover Eq. (35). For this potential, 
the form factor V^,k is given by a rather complicated expres- 
sion, which will not be reproduced here. Important features of 
V^,k are the presence of minima and a decrease with increasing 
asymptotic momentum. This decrease dampens the oscilla- 
tions, such that Vpk, in comparison to the <5-shell form factor, 
decays much more rapidly for large p. This effect becomes 
more pronounced as the width of the potential increases. 

In Fig. |8|(a), the contribution of the two shortest trajecto- 
ries to the ATI spectra is displayed for the Gaussian potential 
d4l|), in comparison to the <5-shell potential. We considered 
the zero-range-potential form factor Vko [Eq. (p2|)]. As in 
the previous figure, there exist additional oscillations, which 
come from V p k. In Fig. |](b), this is clearly shown, for the 
contributions from the second shortest trajectory. For small 
width, as expected, the <5-shell oscillation pattern is practically 
recovered. For the parameter range considered in the figure, 
this holds for a < 0.5. Major differences are present only for 
a > 1.5. As the width gets larger, there is a displacement in 
the minima of the form factor and a suppression of the pho- 



toelectron yield. This suppression is due to the decay of the 
form factor Vp^. Therefore, even when the shell potentials 
are smoothed out, the oscillations survive. Thus, the possibil- 
ity that they are artificially caused by the sharp edges of the 
5-shell potential can be ruled out. 



VII. CONCLUSIONS 

We investigate the influence of the binding potential in 
above-threshold ionization (ATI) for linearly polarized laser 
fields, in terms of quantum orbits, using the uniform approx- 
imation [Eqs. ( ^6| ) and (f28|)]. In this method, the transition 
amplitude is expanded in terms of the collective contribution 
of pairs of orbits rather than individual orbits. No information 
is required beyond the conventional saddle-point approxima- 
tion. This is made possible and, indeed, necessitated by the 
fact that for laser-induced rescattering phenomena the orbits 
naturally come in pairs that nearly coalesce at the classical 
cutoffs, thus rendering the conventional saddle-point approxi- 
mation inapplicable in this energy region. Moreover, the uni- 
form approximation remains valid beyond the classical cutoff 
in the classically forbidden region, where it automatically in- 
corporates the fading out of unphysical saddles beyond the 
cutoff energy. If the two saddles of a pair are sufficiently far 
apart, the standard saddle-point approximation is recovered. 

The fact that the uniform approximation is valid in the 
whole energy range, both away from as well as near the cut- 
offs, allows one to obtain quantitative predictions for ATI 
spectra. Indeed, in this paper this approximation has been 
tested for the zero-range potential against the numerical com- 
putation of the SFA transition amplitudes. The photoelec- 
tron spectra, as well as the angular distributions obtained in 
both ways turned out to be practically identical. With the 
conventional saddle-point approximation, quantitative predic- 
tions are not possible in certain energy regions, which for low 
laser intensity can span the better part of the ATI plateau. 

The excellent quality of the uniform approximation for the 
zero-range potential also suggests that the uniform approxi- 
mation is reliable enough for computing ATI spectra for other 
binding potentials, such as Coulomb, Yukawa, or shell po- 
tentials. Within the framework of this paper, the influence of 
the binding potential is contained in two form factors, which 
either characterize the transition from the ground state to an 
intermediate momentum state, or the transition from the inter- 
mediate state to an asymptotic momentum state. Throughout 
the paper, these form factors are called Vko and V p k, respec- 
tively. 

As a first application, we investigated the role of 
the Coulomb tail by computing photoelectron spectra for 
Coulomb and Yukawa potentials. As a main feature, we 
observe a suppression of the photoelectron yield for the 
ATI plateau, in comparison to the zero-range case, for both 
Yukawa and Coulomb cases. This is due to the functional 
forms of Vpk, which are inversely proportional to the photo- 
electron momentum. Additionally, for the Coulomb poten- 
tial this form factor causes an increase in the low-energy ATI 
peaks. These results are in agreement with the fully numeri- 
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cal solution of the time-dependent Schrodinger equation [|3z 
Furthermore, for the Yukawa potentials, we observed an in- 
crease in the yield for decreasin g sc reening parameter. Similar 
features have been obtained in [22|, from the numerical solu- 
tion of the strong-field approximation transition amplitudes. 

Another class of potentials that we investigated are shell 
potentials, which are commonly used as an approximation for 
clusters. In comparison to the zero-range case, the photoelec- 
tron spectra computed for such potentials exhibit additional 
structure, which comes from the oscillating form of V p k. This 
is an extreme case of how the form factor V p k influences the 
photoelectron yield. Such oscillations are also present when 
the potentials are smoothed out, and therefore are not an arti- 
fact of the shell models. 

An alternative for performing such investigations is the nu- 
merical solution of the three-dimensional Schrodinger equa- 
tion. This would require considerable numerical effort, and, 
for elliptical polarization, it would take one close to the 
limit of today's computational resources. Another possibility 
would be the numerical solution of the strong-field approxi- 
mation amplitudes (jlj) and (Q). From the numerical viewpoint, 
this is not an easy task either, since one must deal with mul- 
tiple integrals of highly oscillating functions. Thus, the uni- 
form approximation considerably simplifies the computations 
involved. Furthermore, using this approximation, one is able 
to gain additional physical insight into the interference pro- 



cesses between the quantum orbits, and how such processes 
are affected by the binding potential. 

Summarizing, the uniform approximation is a very pow- 
erful method for investigating laser-assisted rescattering pro- 
cesses, being applicable in all energy regions of the spec- 
tra. This approximation allows one to compute photoelec- 
tron spectra for binding potentials other than the zero-range 
with minimal numerical effort. Application of the methods 
developed in this paper to other high-intensity laser-induced 
or laser-assisted phenomena, such as non-sequential double 
ionization, or to elliptically polarized fields is, in principle, 
straightforward. 
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